Scintillation of liquid neon from electronic and nuclear recoils 
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We have measured the time dependence of scintillation light from electronic and nuclear recoils in 
liquid neon, finding a slow time constant of 15.4 ±0.2 /is. Pulse shape discrimination is investigated 
as a means of identifying event type in liquid neon. Finally, the nuclear recoil scintillation efficiency 
is measured to be 0.26 ± 0.03 for 387 keV nuclear recoils. 
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INTRODUCTION 

Liquid neon has been suggested as a low-background 
scintillation medium for the CLEAN (Cryogenic Low En- 
ergy Astrophysics with Noble gases) experiment, which 
will have simultaneous sensitivity to both low energy so- 
lar neutrinos and WIMPs % El El- The CLEAN de- 
tector will use approximately 100 tonnes of liquid neon 
surrounded by photomultiplier tubes (PMTs) in a spheri- 
cal geometry. Scintillation light resulting from both elec- 
tronic and nuclear recoils will be detected by these PMTs. 
Liquid neon is a good candidate for low background ex- 
periments because it has no long lived radioactive iso- 
topes, is easily purified, and is dense enough to allow for 
self-shielding, effectively eliminating backgrounds in the 
center of the detector. The solar neutrino signal will con- 
sist of electronic recoils, while a WIMP signal will consist 
of nuclear recoils. The physics goals of CLEAN require 
a method of discriminating between these two types of 
events. 

It was pointed out in Ref. and @ that electronic 
and nuclear recoil events in liquid neon might be distin- 
guished through the use of pulse shape discrimination 
(PSD). Interactions of ionizing radiation with liquid no- 
ble gases produce excimers that can exist in either a sin- 
glet or triplet state. The lifetime of the triplet state in liq- 
uid neon has been previously measured to be 3.9 /is [16|. 
The singlet lifetime is expected to be much shorter, on 
the order of a few nanoseconds, based on measurements 
of other liquid noble gases [1(3] • 

The slow scintillation light emitted by triplet molecules 
is suppressed in intensity by destructive triplet interac- 
tions, primarily Penning ionization and electron-triplet 
spin exchange; these reactions are stronger for higher ex- 
citation densities such as those produced by nuclear re- 
coils. For a given event, measurement of the fraction of 
the total scintillation signal that is in the fast compo- 
nent can be used to determine which type of event has 
occurred. 

PSD has also been studied in detail for liquid xenon 
based experiments P, 0] and can be used to suppress 
gamma ray backgrounds. PSD in liquid helium has also 
been studied in order to separate electronic recoil events 
from 3 He(n,p) 3 H events in the search for the neutron elec- 
tric dipole moment In an extension of the CLEAN 



concept, PSD has been also suggested to counter 39 Ar 
background in single phase liquid argon dark matter de- 
tector s |4| . Previous measurements of liquid and solid 
neon [la ] have compared beta induced scintillation to 5 
MeV alpha induced scintillation, but no studies exist of 
PSD for the nuclear recoils characteristic of a WIMP sig- 
nal. 

Because the WIMP-nucleon scattering rate is expected 
to decay exponentially with nuclear recoil energy, the 
sensitivity of experiments to WIMPs depends crucially 
on energy threshold, which in turn depends both on the 
scintillation yield for nuclear recoils and the light detec- 
tion efficiency of the apparatus. The scintillation yield 
from nuclear recoils relative to that from electronic re- 
coils has been measured recently by several groups for 
liquid xenon [l|, 0, S @|- In contrast, the scintillation 
yield for nuclear recoils in liquid neon has not been mea- 
sured. 

In this paper we describe measurements of the pulse 
shape discrimination efficiency, nuclear recoil scintillation 
efficiency, and time dependence of liquid neon scintilla- 
tion. 



EXPERIMENTAL DETAILS 
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FIG. 1: Schematic representation of the scintillation cell. 

The data presented here were obtained using a 200 ml 
scintillation vessel viewed by two 75 mm diameter photo- 
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multiplier tubes (PMTs) [l8j]. A schematic of the setup 
can be seen in Fig. Q] 

The scintillation vessel is a stainless steel body at- 
tached to two titanium flanges into which 50 mm di- 
ameter sapphire windows are bonded with Stycast 2850 
epoxy. Titanium and sapphire were chosen due to their 
similar thermal expansion coefficients, allowing a leak- 
tight seal when filled with liquid neon. Because liquid 
neon scintillates in the ultraviolet (sa 77 nm) [T^, we 
place inside the vessel secondary windows coated with a 
thin film of tetraphenyl butadiene (TPB 0). The TPB 
shifts the wavelength of the ultraviolet light to approxi- 
mately 440 nm so that it may pass through the windows 
and be detected by the PMTs. The top window is coated 
with 0.29 mg/cm 2 of TPB and the bottom window with 
0.18 mg/cm 2 . The inside wall of the vessel is lined with 
a PTFE sheet that is coated on the inner surface with 
0.16 mg/cm 2 of TPB. 

The scintillation vessel is housed inside a vacuum de- 
war for thermal insulation, and liquid neon is introduced 
though a tube on the side of the scintillation vessel. 
The neon is liquefied from purified gas in a copper can 
mounted to the end of a pulse-tube refrigerator [19(. All 
components that come into contact with the gas or liquid 
are baked to at least 60 C for 12 hours, and the ultra-high 
purity neon gas (99.999%) is passed through a heated gas 
purification getter [20( before it is liquefied. 

A stainless steel tube leading out of the scintillation 
vessel leads to a charcoal trap connected in line. Pu- 
rification of neon with charcoal using this method has 
been previously discussed in Ref. Q. After the cell is 
filled with liquid, a small heater below the trap is ener- 
gized, forcing vaporized neon through the charcoal and 
back into the liquefaction can. By monitoring the heat 
load on the refrigerator we calculate that a continuous 
flow rate of approximately 0.5 gas litres per minute was 
obtained over the course of the experiment. 

Two sample oscilloscope traces from scintillation 
events in liquid neon can be seen in Fig. [51 The sig- 
nal from each PMT is divided three ways with a linear 
fan-out. One copy goes to the triggering electronics and 
two copies go directly to an oscilloscope where they are 
digitized at two different gain settings. We use this tech- 
nique to increase the effective dynamic range of the 8-bit 
oscilloscope digitization. The triggering electronics re- 
quire a coincidence of at least one half of a photoelectron 
(pe) within 100 ns in each PMT. 

All data analysis is performed in software. The first 
analysis step is to take collections of runs, typically two 
hours long, and find the value of the single photoelectron 
signal. For all types of excitations in the neon, there is 
a significant triplet component, and since this is spread 
out over many /lis, a large number of single photo-electron 
events can be analyzed to obtain the gain of the PMTs. 
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FIG. 2: Examples of two events that are identified as elec- 
tronic and nuclear recoils using pulse shape discrimination. 
The discrimination method used is described in the Experi- 
mental Results section. 



Sodium-22 Calibration 




1500 



Electronic recoil signal (pe) 

FIG. 3: Plot of 22 Na calibration data with fit to Eq. JTJ, 
indicating a detector signal yield of 0.93 ± 0.07 photoelectrons 
per keV. 

We use a 10 /LtCi sealed 22 Na source to determine the 
absolute signal yield for electronic recoils. This source is 
placed on the outside of the vacuum enclosure to provide 
a source of high energy gamma rays. Because the gam- 
mas from this source are produced in the e + + e~ -^7 + 7 
reaction, we can significantly reduce backgrounds due to 
other radioactive decays, as well as cosmic rays, by re- 
quiring a coincidence between the scintillation vessel and 
an external organic scintillation detector. The major- 
ity of the 511 keV gamma rays that interact in the neon 
Compton scatter, providing a wide range of energy depo- 
sitions over the region of experimental interest. We also 
use these 'tagged' events for measuring the electronic re- 
coil time dependence. 

Data are analyzed from the two PMTs with the only 
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requirement that the integrated signal from one PMT be 
within a factor of four of the other. We fit the calibration 
data to a function that is comprised of two Gaussian- 
broadened step- functions, 

N = All-erf ( — 

where E 1 = x/F is the energy, Y is the yield, -E511 is the 
value of the 511 keV Compton edge (341 keV) and -B1274 
is the value of the 1274 keV Compton edge (1062 keV). 
<7i and (T2 (= &i v Eyzji) are the widths, and A and 
B are the relative sizes of the steps. The calibration data 
are shown in Fig. [3] along with a four parameter fit for 
<7i, A, B, and Y. From this fit, we have determined that 
the signal yield for electronic recoils in this detector is 
0.93 ± 0.07 pe/keV at 511 keV. 

We measured the stability of the system by taking indi- 
vidual 22 Na runs over the course of four days. We found 
that the signal yield remained constant to within 7%. 

Nuclear recoils 



EXPERIMENTAL RESULTS 
Pulse shape discrimination 

To discriminate between electronic recoil and nuclear 
recoil events, we calculate f prompt, the fraction of the 
signal that arrives in the first 100 ns. 




prompt 



FIG. 5: Plot of the fraction of light in the prompt component 
(fprom.pt) for electronic recoils (Xs) and electronic plus nuclear 
recoils (circles). Events are selected to be between 80 and 120 
photoelectrons for this plot. 




Neon Cell D-D Generator 



FIG. 4: Schematic of the neutron scattering setup. 



Figure [5] is a histogram of f prompt for 1434 events iden- 
tified as electronic recoils and 573 events identified as 
nuclear recoils. All events are selected to be between 80 
and 120 photoelectrons. As one can see, there is a natural 
division at about f pr0 mpt = 0.25. Using a 50% nuclear 
recoil acceptance (f pr ompt > 0.35) we conclude that the 
discrimination efficiency exceeds one part in 1434. 

There is a small population of events in the nuclear 
recoil set that have a small f prompt- The neutron gener- 
ator produces a large number of gamma rays, and some 
of those are misidentified by the time of flight cuts. 



To investigate the detector response to nuclear recoils, 
we use a portable deuterium-deuterium neutron genera- 
tor [2l[ as a source and an organic scintillator as a sec- 
ondary detector. The experimental setup can be seen 
schematically in Fig. |4] The nuclear recoil energy, E rec , 
may be determined from the angle of the scattered neu- 
tron through simple kinematics, 



E r 



2E;, 



[l + Af-cos 2 (6>) 



(1 + Mf 
cos (6)^jM 2 + cos 2 (0) - 1 



(2) 



where Ei„ is the incident neutron energy (2.8 MeV), M is 
the atomic mass (20.18 for neon), and 9 is the scattering 
angle of the outgoing neutron. 



Time dependence 

We estimate the temporal probability density functions 
for nuclear and electronic recoil events in liquid neon 
based on the data from events that yield between 80 and 
120 photoelectrons. Because the d-d neutron generator 
emits significant numbers of gamma rays, we require that 
the fraction of light within the first 100 ns (/prompt) be 
greater than 0.18 before it is considered to be a nuclear 
recoil. There is a clear separation in the neutron and 
gamma populations in Fig. [5] that justifies this cut. Ad- 
ditionally, a 100 ns window in the time of flight spectrum 
is allowed in order to reduce the number of gamma ray 
events in the region of interest. No such cuts are used for 
the tagged 22 Na electronic recoil events. 
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Parameter 


Electronic recoils 


Nuclear recoils 


A 


0.86 ±0.16 


1.33 ±0.13 


B (ns) 


5.0 ±0.2 


3.78 ±0.08 


C 


0.0137 ±0.0004 


0.00230 ± 0.00002 


t (ns) 


-4.8 


±0.9 


n (ns) 


18.6 ± 0.2 


T 2 (jUS) 


15.4 ± 0.2 



TABLE I: Fit parameters to Eq. © for both electronic and 
nuclear recoil events between 80 and 120 photoelectrons. 



Nuclear recoil scintillation efficiency 



FIG. 6: Plot of the probability density function as a function 
of time for both electronic recoils (top curve) and nuclear 
recoils (bottom curve). Solid lines are fits described in the 
text. 



For each event class, the trigger time is defined as the 
time that the voltage rises above 20% of its maximum 
value. These aligned pulses are averaged and then nor- 
malized in order to more easily compare them. Both 
gamma and neutron data are simultaneously fit using a 
least squares method to the following function, 

Pm = Ae-(*- { °)/ T1 
B 



t-to 



(3) 



while allowing A, B, and C to vary between the two 
different types of excitations. 

This function was constructed as an extension of the 
two exponential fit that one needs as a minimal descrip- 
tion of the singlet and triplet decay constants for the neon 
molecules. The quality of fit to the neon data is statis- 
tically better with this hybrid function than it is with 
either a two or three term exponential fit function. The 
motivation for adding the 1/t term is based on observa- 
tions of radiative decay in liquid helium scintillations (l3j 
where the non-exponential component is explained as a 
two-body reaction effect. This is a first order approxima- 
tion to the scintillation time dependence for triplet-triplet 
interactions calculated by King et al. ll| . The resulting 



fit parameters are listed in Table I. We find the slow time 
constant to be significantly longer than that measured by 
Michniak et al. 16] (15.4 ± 0.2 (is vs. 3.9 ± 0.5 (is). We 
attribute this discrepancy to our continuous purification 
procedure which reduces impurities that may quench the 
triplet molecules. 

The singlet decay time can not be directly calculated 
from Ti as the speed of the prompt pulse is also deter- 
mined by the PMT timing and TPB wavelength shifter 
lifetime. 
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FIG. 7: Plot of the light production for 387 keV nuclear recoils 
with Gaussian fit. 

The relative nuclear recoil scintillation efficiency is the 
ratio of light produced by nuclear recoils compared to 
that produced by electronic recoils of the same energy. 
For simplicity we assume that the electronic recoil scin- 
tillation efficiency is constant over our region of interest. 

We choose a scattering angle of 120° ± 3° which 
corresponds to a nuclear recoil energy deposition of 
387 ± 11 keV. We exclude events from the raw data 
set that lie outside of a 60 ns time of flight window. The 
resulting data are plotted in Fig. [7] 

The Gaussian best fit parameters yield a peak at 93 ± 2 
photoelectrons and a width of 33 ± 3 photoelectrons. 
Therefore, the nuclear recoil signal yield is 0.24 ± 0.01 
photoelectrons per keV. Using 0.93 ± 0.07 pe/keV for 
511 keV electronic recoils, the relative nuclear recoil scin- 
tillation efficiency is 0.26 ± 0.03 at 387 ± 11 keV. 



CONCLUSION 

We have measured the time dependence of electronic 
and nuclear recoil scintillation events in liquid neon. We 
find that the slow time constant is 15.4 ± 0.2 /is and find 
evidence for a 1/t term through least squares fitting. We 
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have also found that in our apparatus the signal yield 
for 511 keV electronic recoils is 0.93 ± 0.07 pe/keV and 
for 387 keV nuclear recoils is 0.24 ± 0.01 pe/keV. This 
implies a relative scintillation efficiency of 0.26 ± 0.03 at 
this energy. For a signal yield between 80 and 120 photo- 
electrons we find the pulse shape discrimination between 
nuclear and electronic recoils to be better than 1 part in 
1434 with 50% nuclear recoil acceptance. 
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